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Abstract 

Taking advantage of the electron Bernstein waves heating (EBWH) 
system of the TJ-II stellarator, an electron Bernstein emission (EBE) 
diagnostic was installed. Its purpose is to investigate the B-X-0 radi- 
ation properties in the zone where optimum theoretical EBW coupling 
is predicted. An internal movable mirror shared by both systems al- 
lows to collect the EBE radiation along the same line of sight that is 
used for EBW heating. The theoretical EBE has been calculated for 
different orientations of the internal mirror using the TRUBA code 
as ray tracer. A comparison with experimental data obtained in NBI 
discharges is carried out. The results provide a valuable information 
regarding the experimental O-X mode conversion window expected 
in the EBW heating experiments. Furthermore, the characterization 
of the radiation polarization shows evidence of the underlying B-X-0 
conversion process. 



1 Introduction 

Plasma heating by Electron Bernstein waves [1] has been successfully demon- 
strated in several magnetic confinement devices [2, 3, 4, 5]. The main advan- 
tage of this technique is that it offers the possibility to heat plasmas above 
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the density limits of the standard electromagnetic modes. For a complete 
review of the EBW topic see [6] and references therein. Moreover, the emit- 
ted Bernstein waves from the over-dense plasma core can be helpful for the 
plasma diagnosis in high density scenarios. See, for instance, ref. [7], where 
EBE is used as an electron temperature diagnostic and ref. [8] where it is 
employed to measure the magnetic field. 

The TJ-II stellarator [9] is a medium size heliac with major radius Rq = 
1.5 m, minor radius a ~ 0.2 m and magnetic field strength on axis around 
0.95 T. In this device, the use of the available ECR heating sources -two gy- 
rotrons delivering 300 kW power each and injecting X-mode polarized waves 
at second harmonic (/ = 53.2 GHz)- is limited by the cut-off density for 
that frequency and that mode, i.e. n e m 1.7 x 10 19 m~ 3 . Two neutral beam 
injectors allow to heat the plasma above this limit and therefore, additional 
ECRH heating in NBI plasmas becomes possible only through electron Bern- 
stein waves. The heating scheme using Bernstein waves for TJ-II is based on 
the O-X-B double mode conversion using low field side launching [10]. 

As it is well-known, the efficiency of the O-X mode conversion depends 
strongly on the launching direction and the polarization of the incident wave. 
Thus, in order to find an optimum experimental heating efficiency, a steer- 
able mirror was installed inside the vacuum vessel, allowing a horizontal and 
vertical scan of the ECRH beam launching direction around the direction for 
which maximum O-X mode conversion efficiency is expected. Note that a 
perfect optimization procedure should also include the displacement of the 
mirror center, since for each penetration location of the plasma wave a dif- 
ferent well defined injection direction is required. This is not possible with 
the existing design. 

In the present setup, the mirror can also be positioned so that the plasma 
radiation coming along the EBWH launching direction is redirected towards 
a radiometer antenna that measures thermal emission at 28 GHz. With this 
setup, measuring and heating can not be performed at the same time. Since 
both the EBWH and the ECE/EBE systems share the same line of sight, 
the analysis of electron cyclotron radiation in plasmas above the O mode 
cut-off density for 28 GHz (n e ^ 1.0 x 10 19 m 3 ) may provide a very valuable 
information to determine the optimum direction for O-X-B heating. As a 
matter of fact, because of the low cut-off densities of the 28 GHz electro- 
magnetic modes, any radiation detected in plasmas above this density must 
come from an inverse B-X-0 mode conversion unless the contribution of the 
peripheral underdense plasma to the emission is high enough to mask the 
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B-X-0 component. This point will be addressed in the discussion. 

In this work, we present the theoretical characterization of the EBE re- 
sulting from B-X-0 conversion, for different orientations of the launching 
mirror and different rotation angles of the radiometer, in the vicinity of the 
position that provides the maximum theoretical O-X conversion efficiency. 
The numerical results are obtained with the ray tracing code TRUBA and 
the well-known concepts in traditional ECE calculations, i.e. energy balance 
between emission and absorption as well as black body approximation. The 
theoretical results are compared with the experimental ones, obtained under 
ECRH + NBI discharges. 

The article is outlined as follows: section 2 focuses on the description of 
the experimental setup; section 3 deals with the numerical EBE calculations; 
section 4 is devoted to the experimental results and their comparison with 
the predictions presented in section 3; and finally, section 5 summarizes and 
discusses the main results. The appendix describes the calculation of the 
expected power in the different polarizations measured by the radiometer 
when it is rotated around its symmetry axis. 

2 EBWH launching system and EBE diag- 
nostic 

The EBWH launching system consists of a corrugated waveguide that trans- 
mits up to 300 kW of microwave power into the vacuum vessel, where a 
steerable ellipsoidal mirror (170 mm xl90 mm) focuses the beam and defines 
the launching direction. Prior to coupling the radiation to the waveguide, 
two elliptical focusing mirrors and two A/4 and A/8 plane polarizers provide 
a Gaussian beam with the properly polarized electromagnetic field. For the 
details of the EBW heating system see refs. [11, 12]. The position of the 
movable mirror is determined by two orientation angles along the horizontal 
(or toroidal) and vertical (or poloidal) directions, and ag respectively. If 
the launching direction is given by these angles, the radiation coming from 
this same direction can be redirected to the 28 GHz heterodyne radiometer 
using a different pair of angles a.'^ and a' e and a flat mirror attached to a 
corrugated waveguide, similar to the one used in the launching system (see 
figure 1). After its passing through a glass focusing lens located outside the 
vacuum chamber, the radiation is detected by a quad-ridged dual-polarized 
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Figure 1: EBWH launching system and EBE diagnostic. Different mirror posi- 
tioning angles are used for EBW heating (a^, atg) and for EBE detection (a',, a' e ). 
Radiated power suffers two consecutive reflections and is coupled to the detection 
waveguide. 



microwave horn. The lens modifies the beam pattern at the output of the 
waveguide in order to obtain an efficient coupling to the horn antenna. For 
further details of the EBE diagnostic see [13, 14]. 

Once the 28 GHz radiation reaches the horn, it is split into two sepa- 
rate components according to two orthogonal detection directions. Properly 
calibrated, the sum of both signals, 7 EBE i an d / E be2, corresponds to the to- 
tal radiative temperature T rad . Accordingly to the properties of an obliquely 
propagating O mode, an elliptical polarization is expected. Therefore, in 
order to perform an experimental measurement of this polarization, the ori- 
entation of the detection directions may be changed from shot to shot by 
rotating the diagnostic an angle £ around its symmetry axis. Actually, only 
the amplitude of the field components can be measured to some degree of ac- 
curacy and no information about the interacting phase between components 
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Figure 2: (a) The rays are launched perpendicularly to the wave front surface at 
the plasma edge. 



and thus about the rotation sense of the polarization can be extracted. 

Using purely geometrical calculations with the vacuum field configuration 
it can be seen that for ( = 10°, channel 2 (EBE2) is aligned along the 
theoretical major axis of the expected O mode polarization ellipse. Channel 
1 is aligned along the minor axis of the ellipse. See the Appendix for a more 
detailed description of the diagnostic rotation and the expected polarization. 

3 EBE numerical simulations 

The calculation of the B-X-0 radiation has been performed with ray tracing 
simulations using the TRUBA code, whose detailed description can be found 
in refs [15] and [12]. 

From now on, the internal mirror angles that we set to collect the 28 GHz 
radiation will be referred as and ag, dropping the primes. 

To simulate the emission, we need first the ray trajectories obtained by 
tracing a beam launched from the corrugated detection antenna towards the 
plasma (all the simulations have been done launching 121 rays, and for the 
standard equilibrium configuration, provided by VMEC-based libraries [16]. 
The maximum averaged /3 never exceeds 1%, thus no significant modification 
in the equilibrium is produced, since TJ-II is designed to present a very small 
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Shafranov shift. The ellipsoidal mirror has been designed to focus the beam 
generated by the launching antenna. The incidence angles and the distance 
to the detection waveguide for EBE differ from the ones used in the EBW 
heating configuration and therefore, the mirror is no longer optimized for 
the EBE measurements. The changes that this lack of optimization may 
produce are neglected here. Moreover, note that the ray tracing method has 
a fundamental limitation when it is applied to the O-X conversion process 
since it cannot take into account the beam spectrum in a fully self-consistent 
way. Therefore, the total O-X conversion efficiency obtained from adding 
the different contribution of each ray -properly weighted with the Gaussian 
beam profile- differs in general from the O-X conversion efficiency of the real 
electromagnetic field of the beam [17]. We will discuss this point together 
with the interpretation of the emission data when we come to the summary. 

Figure 3 shows the result of a typical ray tracing simulation. The electron 
density and a temperature profiles used are fitted to the Thomson Scattering 
(TS) ones obtained in NBI-heated TJ-II plasmas. To calculate the electron 
Bernstein emission, the radiative transfer equation [18] is solved along the tra- 
jectory of the O-X transmitted rays. The solution to this equation provides 
the emission intensity per unit frequency, which is given by the well-known 
expression 



where T e is expressed in energy units, c is the speed of light, oj = 2irf is the 
wave angular frequency and dr = — a^da is the optical thickness, with a u the 
absorption coefficient and a the ray coordinate. For each transmitted ray, the 
integral in eq. 1 is performed along the ray path up to the inner boundary 
of the conversion layer, i.e. the SX-mode cut-off, in order to obtain the 
emission intensity Ii prior to the single ray X-0 tunneling process. Taking 
into account the conversion efficiency for each particular ray, r]i, and the 
Gaussian beam profile weighting factor we may write the total emission 
intensity as 



where i is the ray index and n c is the total number of converted rays. 




(1) 




(2) 
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Figure 3: Ray tracing with 121 rays (only 21 rays are represented). The beam 
transmission efficiency is around 90%. The black solid lines correspond to the 
central trajectories of the launched, transmitted and reflected rays (see legend). 



3.1 Dependence on mirror position, density and tem- 
perature 

In order to compare with the result of the experiments presented in section 4, 
the EBE calculations have been carried out for different lines of sight. The ex- 
perimental TS profiles have been fitted to the functions n e {p) = no(l— p 2ni ) n ' 2 
and T e (p) = T (l — p 2 * 1 )* 2 , where the set of parameters {^1,^2,^1,^2} = 
{1.6,2.7,2.0, 1.7} is found to match best the whole set of experimental pro- 
files. A wide range of the no and To values have been considered. Figure 
4 shows some of the experimental and analytical electron density (a) and 
temperature (b) profiles. The dependence of the electron Bernstein emission 
on the mirror positioning angles is plotted in the figures 5 for a fixed central 
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Figure 4: Experimental electron density (a) and temperature (b) Thomson scat- 
tering profiles obtained in NBI heated plasmas. The numerical calculations have 
been carried out using a fit to the average shape of these profiles. 



density and four different values of central temperature T . As expected, the 
emission intensity increases with the electron temperature. In all cases, the 
maximum emission is located around ~ 31° and ag ~ 21.5°. 

It is observed that the emission finds a stronger dependence on ag than 
on a<f,. This is due to the more pronounced curvature of the plasma along the 
poloidal direction. On the one hand, this makes the incidence direction at 
the conversion layer deviate more sensitively from its optimum ag value than 
from the a^ one, which affects to the conversion efficiency. And on the other 
hand, it changes appreciably the optical path of the rays and consequently 
the radiation reached at the O cut-off layer previous to conversion. 

Regarding the dependence of EBE on central density, note that, as n 
increases, the O-mode cut-off layer moves radially outwards. This modi- 
fies the O-X conversion efficiency [19] through its dependence on the values 
at the O mode cut-off of the characteristic density gradient scale length 
(L n = n e /|Vn e |), Nn and N±, and the magnetic field strength B. The final 
combination of these effects is illustrated in figure 6, where, independently of 
the central temperature To, a maximum intensity is observed always around 
n = 3 x 10 19 m~ 3 . 
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Figure 5: Dependence of the 28 GHz B-X-0 mode converted radiation intensity 
for (a) T = 0.1, (b) T = 0.3, (c) T = 0.5 and (d) T = 0.7 KeV. In all the 
cases, a central density given by uq = 3.5 x 10 19 m~ 3 has been used. Each map 
represents the result obtained with 256 launched beams and 121 rays per beam. 



3.2 Dependence on the rotation angle £ 

By rotating the diagnostic we should be able to determine if there is a dom- 
inant polarization in the incoming radiation and if this coincides with the 
expected one. For the ideal case of a perfect elliptically polarized wave com- 
ing from the B-X-0 emission, the power dissipated along each detection 
direction when the diagnostic rotation angle is set to some arbitrary ( is 
given by eqs. 8 of the Appendix. In order to consider the presence of an 
unknown amount of non-polarized radiation we define 
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Figure 6: Dependence of the B-X-0 radiation intensity on no and To for the line 
(a^ag) = (31°, 21.5°). 



where P& = fdP2(C = 10°) is taken as a fraction of the maximum amplitude 
measured in channel 2. A minimum in 11(C) is expected for C — 10°, when 
channel 2 is directed along the major axis of the polarization ellipse. The 
dependence of 1Z on the rotation angle C of the horn antenna is represented 
in figure 7 for different values of fa- As it is clear, the greater the amount of 
unpolarized radiation, the flatter the 1Z(C) profile becomes. 



4 Experimental results 

Figure 8 shows the time evolution of the main plasma parameters observed 
during a typical ECRH+NBI shot. Plasma heating is first performed by 
applying 500 kW of second harmonic off-axis ECRH power. Then, NBI 
heating (P NB i ~ 1 MW) is turned on and densities in the range 2 — 3.5 x 10 19 
m -3 are achieved. ECRH power is turned off before the second harmonic cut- 
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Figure 7: Dependence on the power ratio TZ{Q on the EBE diagnostic rotation 
angle. Different non-polarized background radiation fractions are considered 



off density for 53.2 GHz (n e « 1.7 x 10 19 m -3 ) is attained. This provokes the 
drop in the ECE signals, which still show a low plasma temperature until the 
NBI-induced electron density growth prevents the detection by ECE means. 
The increase in the electron temperature and density during the NBI phase 
is accompanied with the subsequent raise of radiation observed by the Soft 
X ray (SXR) detector. EBE and ECE radiation after the ECRH turn-off are 
compared in figure 9. The ECE signal decreases faster because this diagnostic 
is looking to the plasma along a perpendicular line of sight whereas the EBE 
diagnostic is looking with a highly oblique view. Thus, the 28 GHz emission 
comes from suprathermal electrons with a lower collisionality than that of the 
thermal electrons responsible of the emission detected by the ECE diagnostic. 

4.1 Experimental and ag scans 

Figures 10 and 11 show the measured radiative temperature at the TS ac- 
quisition time for different positioning angles of the internal mirror. The 
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Figure 8: Typical time evolution of some relevant quantities obtained in 
ECRH+NBI heated plasmas. 
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Figure 9: Zoom view of the central ECE and EBE2 signals represented in figure 
8. 



first of them corresponds to a scan in a^ for ag = 21.5°, while the second 
one represents, for a^ = 31°, the scan in ag. In both figures T rad is referred 
to the left axis (blue dots). In order to remove its strong dependence on 
the temperature profile, which was not constant between shots, these values 
are normalized to the central temperature value {Tq S ) provided by the TS 
diagnostic (red dots, right axis). 

The scanned positions cover a wider area than the one considered in the 
theoretical calculations. The result presented in figure 10 is consistent with 
the theoretical predictions. The measured intensity show an almost flat pro- 
file weakly dependent on a^. Regarding the scan shown in figure 11a stronger 
variation of the collected 28 GHz intensity on ag is clearly manifested, as was 
also expected. The value for which a maximum emission value is obtained is 
ag ~ 19.5°. This can be interpreted as an acceptable result compared with 
the theoretical prediction {ag ~ 21.5°), considering that the range of move- 
ment of the mirror along the poloidal direction is 30° approximately. On 
the other hand this disagreement of approximately 2° represents a noticeable 
deviation compared with the width of the EBE window along the poloidal di- 
rection, which has 4° (see figures 5). Furthermore, if all the radiation received 
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Figure 10: Dependence of T ra d on for a# = 21.5°. The shadowed area represent 
the interval considered in the theoretical emission plots. 



came only from B-X-O-converted one, the intensity should drop drastically 
out of the 4° of width of the viewing window. Since this is not what happens, 
it can be concluded that the amount of 28 GHz disperse radiation collected 
by the mirror is far from being negligible. Such a contribution of disperse 
radiation is probably coming from EBE emitted by the over dense plasma 
core all around the torus. This radiation leaves the plasma via B-X-0 or 
B-X conversion and reaches the antenna after reflecting on the walls, which 
also provoke randomly changes on its polarization. 

4.2 Experimental determination of polarization 

A scan of the rotation angle of the diagnostic has been carried out in order to 
verify that the radiation reaching the mirror contains a direct contribution 
of O-mode polarized radiation. If this contribution exists, we should observe 
a similar trend in the experimental and the theoretical ratio 7Z{() (see fig. 
7). The polarization has been measured along the line of sight of the mirror 
for which the measured EBE intensity is maximum along the poloidal angle 
(see fig. 11). This line of sight corresponds to the values olq = 19.5° and 
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Figure 11: Dependence of T ra( j on a# for = 31°. 



= 31°. The rotation angle ( was varied from —10° to 90°. Figure 12 
shows the experimental ratio 11(C) = ^ebei(C)/^ebe2(C) measured at the TS 
acquisition time and the calculated one when fa = 1, that is, the value for 
which the best agreement is observed. Taking also into account the result 
shown in fig. 11, it can be concluded that the background unpolarized field 
contributing to the measured intensity is of the same order of magnitude of 
the O-mode radiation reaching the detector. 



5 Summary and discussion 

In this work a comprehensive study of the 28 GHz radiation rising from B-X- 
O conversion in over-dense TJ-II NBI-heated plasmas has been undertaken 
theoretically and experimentally. The numerical results have been obtained 
with the TRUBA ray tracing code. As expected from an X-0 conversion 
process, highly dependent on the wave propagation direction, a strong vari- 
ation of the intensity has been found as the viewing angle along the poloidal 
direction is modified. It has been observed that the experimental maximum 
of the emitted radiation was shifted 1.5 — 2.0 degrees with respect to the 
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Figure 12: Experimental ratio 1Z(C) = / B bei/^ebe2 compared to the theoretical 
estimation, for = 1. A significant fraction of disperse radiation is found. 

numerical estimations, which is not negligible compared with the approxi- 
mately 4 degrees of width of the EBE window along the poloidal direction. 
As a task for a forthcoming work, a fine scan on this angle should be per- 
formed in order to resolve that maximum. Regarding the maximum value 
of J EBE along the toroidal direction, the observation is in agreement with the 
theoretical prediction, since both theory and experiment shows an almost 
flat profile / BB e(o!</p)- From these results, an optimum experimental direction 
for forthcoming EBW heating experiments is extracted. 

The polarization of the collected radiation has been measured in order 
to verify that the received radiation comes mainly from a B-X-0 conversion 
process. The experiment has confirmed the existence of the expected O-mode 
component, importantly screened by disperse radiation. 

Further discussion can be made on the possibility that this O-mode could 
be emitted from the under-dense plasma edge. This point is clarified in 
figs. 13(a) and 13(b), which show the O-mode absorption coefficient for 
different values of Nn in a cylindrical plasma with typical TJ-II magnetic 
field, electron density and temperature profiles. This has been estimated by 
using well-known expressions of the weakly relativistic absorption coefficient 
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Figure 13: (a) Absorption coefficient (solid lines) for oblique O-mode propagation 
at 28 GHz as a function of the effective radius, considering a cylindrical plasma 
with the magnetic field profile similar to the TJ-II one. The electron temperature 
(dashed line) and density (dotted line) profiles like those generated in ECRH- 
heated TJ-II plasmas are also represented. The density is under the O-mode 
cut-off in all the profile, (b) Oblique O-mode absorption coefficient for the over- 
dense case, i.e. profiles from NBI-heated TJ-II plasmas are considered. In both 
cases the different values of Nu considered are in the range between Nu ~ and 
N\\ ~ 0.65. The more peaked curves of the absorption coefficient correspond to 
the lower value of N\\ . 



of electromagnetic modes [20]. On the one hand, figure 13(a) represents the 
case where the density is below the O-mode cut-off, i.e. between 1025 and 
1100 ms in figure 8. On the other hand, figure 13(b) shows the over-dense 
case in the interval between 1100 and 1150 ms. The absorption coefficient in 
the over-dense scenario is three orders of magnitude lower than in the under- 
dense case. Since in this situation the absorption of the O-mode is negligible 
in the edge, so is the emission coefficient. This allows to conclude that the O 
mode polarized radiation measured when in over-dense conditions can only 
be attributed to a B-X-0 process. 

A considerable amount of disperse radiation has been found. Its presence 
can be understood by considering that the emission of EBWs takes place all 
around the over-dense plasma column and reaches the antenna after many 
reflections in the vacuum chamber, which in turn produces a random po- 
larization. Further studies on the characterization of the disperse radiation 
using ray tracing are planned. This task results extremely important due to 
the role of this contribution to the deterioration of the EBE signal. As a 
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matter of fact, the disperse radiation together with the presence of the X-0 
conversion layer in between the emitting source and the radiometer horns, 
represents a handicap in the calibration of the EBE signal by traditional 
means if the electron temperature diagnosis in over-dense regime is aimed. 
Note that X-0 conversion layer act as a kind of k-spectrum filter that can 
not be taken into account in the calibration procedure. It could only be in- 
cluded a posteriori knowing the exact O-X conversion efficiency of the beam 
and assuming a symmetric conversion efficiency, which is not always true. 
Moreover, as it was briefly mentioned in section 3, the calculation of the 
beam O-X conversion efficiency by ray tracing is not accurate, making the 
final determination of an electron temperature a very hard task. 

Finally, the improvement in the spatial resolution by considering a more 
collimated viewing pattern as in standard ECE measurements could help to 
get rid of the disperse radiation. However, this will be accompanied by a 
X-0 conversion efficiency drop due to the widening of the k-spectrum of 
the field associated to a highly collimated beam and, consequently, the EBE 
source for the temperature measurement would be strongly screened by the 
disperse radiation. 
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Appendix: Expected power for an arbitrary 
rotation angle 

For the derivation of an explicit form of the powers Pi(£) and ^(Oj the 
ideal theoretical polarization of the B-X-0 radiation is first defined in the 
main reference system given by the polarization vectors ei and e 2 (see figure 
14). As already anticipated, the azimuth of the O mode polarization ellipse 
in this reference system can be obtained from geometrical calculations by 
finding the projection of the {k, B} plane, defined at the plasma boundary, 
onto the detection plane of the horn, after the two successive reflections both 
in the elliptical and the flat mirror. Furthermore, in order to couple an 
obliquely propagating ordinary mode at the plasma edge the wave must be 
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Figure 14: Theoretical O-mode polarization ellipse. In the figure, the azimuth is 
defined in the non rotated frame (£ = 0). The polarizations vectors of the detector 
when it is rotated an angle e EBE1 (C) and 

Gebe2(C)> are a ^ so represented. 



elliptically polarized with an ellipticity angle 7 given by 

tan 7 =^- — ^, 4 

2 cos 9 

where g 2 = sin 4 9+{A/Y 2 ) cos 2 9, Y = co ce /uj and 9 is the propagation angle at 
the plasma edge. For the standard TJ-II magnetic configuration and taking 
into account the fact that the detector is looking to the emitted wave (the 
rotation sense of the electric field is reversed), we finally obtain the angles 
that define the polarization ellipse in the main reference system, i.e. 

ij Q = 100°, xo = -7 = "35°. (5) 

Taking \Ei\ 2 + \E 2 \ 2 = 1 the relation between these angles and the wave 
electric field is given by [21] 

si = E 2 — E 2 = cos 2xo cos 2ip Q , 

s 2 = 2|_E 1 ||_E 2 | cos5 = cos2x sin2^ , (6) 
S3 = 2 1 Ei 1 1 Ei I sin 5 = sin2x , 

where |-Ei|, I-E2I and 5 = S 2 —Si determine the complex electric field amplitude 
through E = Eiei + E 2 e 2 = \Ei\ exp(i<5i)e! + \E 2 \ exp(i<5 2 )e 2 . Consider now 
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the field components in an arbitrary (^—rotated system defined by the two 
orthogonal directions of the EBE diagnostic: E = E[e EBE1 + P 2 £EBE2-Note 
that when ( « 10°, the vector e EBE2 is directed along the major axis of the 
polarization ellipse. Therefore, the power detected in both channels, for an 
arbitrary angle (, is given by 

Pi(C) = E[E[* = |£x| 2 cos 2 C + \E 2 \ 2 sin 2 C + 2^(E 1 E 2 ) cosCsinC 
P 2 (C) = E' 2 E' 2 * = |E 2 | 2 cos 2 C + \E 1 \ 2 sm 2 (-2$l(E 1 E 2 )cos(sm(. 

Recovering the relations (6) and taking 5\ = 0, Pi and P 2 can be written as 

sin 2 ( + s 2 cos ( sin £ 

(8) 

sin 2 ( — s 2 cos ( sin (. 



l + 2 , . ( 1 - S] 



2 y * v 2 



p 2 (C) - cos ° 

For C = 10°, P 2 is maximum. 
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